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ABSTRACT 

We present observations of CH3CN (12—11) emission at a resolution of 
~2" toward the shell-like ultracompact Hii region G5. 89— 0.39 with the Submil- 
limeter Array. The integrated CH3CN emission reveals dense and hot molecular 
cocoon in the periphery of the Hii region G5. 89— 0.39, with a CH3CN deficient 
region roughly centered at G5. 89— 0.39. By analyzing the CH3CN emission using 
population diagram analysis, we find, for the first time, a decreasing temperature 
structure from 150 to 40 K with the projected distance from Feldt's star, which 
is thought to be responsible for powering the H 11 region. Our results further indi- 
cate that the majority of the heating energy in the observed dense gas is supplied 
by the Feldt's star. From the derived CH3CN column density profile, we conclude 
that the dense gas is not uniformly-distributed but centrally-concentrated, with 
a power-law exponent of 5.5 for r < 8000 AU, and 2.0 for 8000 AU < r < 20000 
AU, where r is the distance to Feldt's star. The estimated large power index 
of 5.5 can be attributed to an enhancement of CH3CN abundance in the close 
vicinity of Feldt's star. 

Subject headings: Hii regions — ISM: individual (G5. 89— 0.39) — ISM: clouds 
— stars: formation 



1. Introduction 

G5. 89— 0.39 (hereafter G5.89, also known as W28 A2) is a shell-like ultracompact (UC) 
Hii region with an angular diameter of ~4", presumably powered by a young 0-type star 
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( IWood &: Churchwel]||l989l ). Recently, nea r-IR observations reported a candidate for power- 
ing the UC Hii region, an 05-V type star (IFeldt et al.ll2003l . hereafter Feldt's star), although 
there is a positional offset of ~1" between the Hii region center and t he near-IR source . 
The off-center location probably results from the migration of Fe ldt's star (IFeldt et al.ll2003[) . 



The estimated distance to G5.89 varies from 1.9 kpc to 3.8 kpc ( jHunter et al 



20081 . and ref- 



erences ther ein). In this letter, we adopt a distance of 2 kpc, which is favored by most recent 



studies (e.g.. lWatson et al.ll2007l : iHunter et al 



mation such as energetic ou tflows (e.g. 



20081). Several signposts of high-mass star for 



SoUins et al 



200 



3 



Watson et al. 2007: Hunter et al. 



20081 ) and maser activities (IFish et al.l l2005l : iKurtz et al.l l2004j ) have been reported toward 
this region. The orientation of the reported outflows in various tracers is notably different 



indicating the existence of multiple young stellar objects in this region (IHunter et al.l 12008 



and reference therein). Indeed, subarcsecond sub-millimeter observations have identified 
at least five dust co ndensations (denoted as SMAl, SMA2, SMA-N, SMA-E and SMA-S) 
( jHunter et al.ll2008l ). With the existence of a young 0-type star still embedded in its na- 
tal cloud core, G5.89 provides an ideal laboratory for studying the physical and chemical 
conditions of molecular gas surrounding a newly formed massive star. 

In this letter, we present methyl cyanide (CH^CN) observations of G5.89 at an angular 
resolution of ~2" with the Submillimeter Arra}t| (SMA). As CH3CN is a symmetric-top 
molecule, its /^-ladder components between adjacent J levels are closely spaced in frequency 
while dramatically different in excitation energies. Therefore, an excitation analysis of the 
CH3CN iT-ladder transitions provides an ideal tool for pro bing (rotational) temp erature 
of dense molecular gas with small calibration uncertainties (ICummins et al.lll983l ). Since 
radiative transitions are not allowed across /^'-ladders, the derived rotational temperature 
is close to the kinetic temperature if th e gas is thermalized. Using population diagram 
(PD) analysis (IGoldsmith fc Langerlll999t ). we derive the distribution of the temperatures as 
well as the filling factor and column density of CH3CN gas toward G5.89. We then discuss 
the implications of the observed temperature, filling factor and CH3CN column density 
structures and evaluate the density profile in this region. A possible variation of CH3CN 
fractional abundance is also illustrated. 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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2. Observations and Data Reduction 

The observations of CH3CN (12—11) in tlie band of 220 GHz were carried out witli 
the SMA on 2008 April 17 and May 2. With seven antennas in the array, the projected 
baselines ranged from about 9 m to 120 m (7 to 90 kA). The phase center was R.A. = 
18'^00™30.32^ (J2000) and decl. = -24°04'00.50" (J2000), and the half-power width of the 
SMA primary beam was ~54". The spectral resolution was 0.41 MHz, corresponding to a 
velocity resol ution of ~0.55 km s~^. The total available double-sideband bandwidth was 



4 GHz. See iHo. Moran. fc Lol (120041 ) for more complete specifications of the SMA. The 
flux calibrators were Uranus (Apr 17) and Titan (May 2), and the bandpass calibrators 
were 3C273 (Apr 17 and May 2) and 3C454.3 (Apr 17). The nearby compact radio sources 
1733-130 {S ~ 2.8 Jy in both Apr 17 and May 2) and 1924-292 {S ~ 4.4 Jy in Apr 17 and 
4.7 Jy in May 2) served as complex gain calibrators. We calibrated the data using the MIR 
software package and made maps using the MIRIAD package. With uniform weighting, the 
synthesized beam size was about 3.1" x 1.8" at P.A. of 55.0°. The rms noise level in a single 
channel of the spectral line images was ~130 mJy beam~^ (equivalently ~0.6 K). 



3. The Spectra and Morphology of the CH3CN Emission 



Figure [T] shows the integrated K = 3 component of the CH3CN (12—11) transition. Fig- 
ur e[2]shows the CH3CN (12—11) spectra toward the five sub-mm dust components identified 
by iHunter et al.l (120081 ) as well as toward Feldt's star. The K components of the CH3CN 
(12—11) transition are detected toward G5.89 up to = 7 with an excitation energy of 
~420 K. Comparing the int egrated CH3CN flux m easured from our SMA observations with 
the single-dish results (e.g.. iPankonin et al.l 1200 ll ). we estimate that approximately 80% of 
the CH3CN emission is recovered by the SMA observations. 

The morphologies of the lower K (i.e., K = 0, 1, and 2) components are similar to 
that of the K = 3 component shown in Figure [TJ A cavity roughly centered at the UC 
Hii region can be discerned. This ca vity is spatially coincident with the dust-deficient 
region reported by iHunter et al.l (120081 ) , indicating a deficiency of both molecules and dust 
grains. With a synthesized beam of 3.1" x 1.8", the molecular cocoon seen in the K = 
3 component is resolved into at least three peaks, here referred to as N-lobe, E-lobe, and 
S-lobe, whose locations agree with the sub-mm dust condensations SMA-N, SMA-E, and 
SMA-S, respectively, as shown in Figure [H The higher K (i.e., K = 5, 6, and 7) components 
are mainly concentrated in N-lobe (i.e., the vicinity of SMA-N as well as SMAl) and become 
much weaker or even undetected toward E-lobe and S-lobe. 
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4. Excitation Analysis of the CH3CN Gas 
4.1. Methodology 

Given the fairly strong detection of CH3CN transitions toward G5.89, estimates of 
gas temperature as a function of position are feasible. We first performed Gaussian line- 
profile fitting to estimate line parameters, including radial velocity, line- width, and integrated 
intensity, of each K component at each pixel. By assuming observed transitions to be 
optically thin and in LTE, we then calculated N'^^, the CH3CN column density of the 
up per level of the transition , from the integrated intensities. For details, see equation (1) 



m 



Thompson fc MacDonaldl (119991 ). Since the inferred gas density of >10^ cm ^ (see 



i s at least a facto r of 10 higher than the critical density at the derived gas temperature 
(jChen et al.ll2006l ). the LTE approximation should be valid. 



In the G5.89 region the lower K components of the CH3CN (12—11) emission are appar- 
ently optically thick, as shown in Figure [2l Under such conditions, a good estimation of gas 



(HoUis 


1982; 


Turner 


1991) 



the PD analysis instead. The difference between the PD analysis and the rotational diagram 
analysis is that the former includes two additional correction factors, i.e., the optical depth 
correction factor CV (= t/(1— e""^)) and the filling factor / {=Qs/^a), where r is the line 
opacity, fig is the source solid angle, and fla is the beam solid angle. Under LTE conditions, 
the rotational temperatures as well as the filling factors and column densities of CH3CN can 
be determined from the following equation in PD analysis: 



Ty-obs 
9n 



Ini 



Q{Tr 



rot J 



kTrni 



+ ln{f) - ln{Cr) 



'V. 



where is the expected upper-state column density of the target molecule incorpo- 
rating the effects of the line opacity (CV) and the beam dilution (/); is the degeneracy of 
upper state; A^'tot is the total column density of the molecule; (^(T'l.ot) is the dimensionless 
rotational partition function; Trot is the rotational temperature; and is the upper energy 
level. 

According to equation (1), for a given upper level, N'^^ can be evaluated from a set 
of Atot, Trot, /, and CV. Since CV is actually a function of Atot and Trot, the independent 
parameters are therefore Atot, T'rot, and / only. We calculated A"*^*^ for the parameter space 
of r,ot = 10-300 K, Atot in the range of 10^3_xoi8 ^^^-2^ y between 0.01 and 1.0. By 
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comparing jSf^^ with N'^^ deduced from observational results, we performed minimization 



X 



(2) 



to determine for each pixel the physical properties, i.e.. Tret, Ntnt, and /, fr om the selected 
parameter space. 6N^^^ is the 1-a error of N'^'^. See IWang et al.l (120091 ) for a detailed 
description of the analysis method. 



4.2. Temperature & Illumination of Molecular Gas 



As shown in Figure [3]a, we have for the first time resolved the temperature structure 
of the neutral dense gas in this region with the PD analysis. The inferred CH3CN gas 
temperatures range from 40 K to 150 K, in agreement with the single-dish results of 56— 75 if 
(IThompson fc MacDonaldlll999l : iPankonin et al.ll200ll : iPurcell et al.ll2006l ) . Among the three 
CH3CN lobes, N-lobe is the hottest, E-lobe next, and S-lobe is the coldest. In particular, 
the CH3CN gas in the close vicinity of Feldt's star has the highest temperature within the 
whole region. 

To further explore the role of Feldt's star in illuminating its surrounding gas, we plot 
the deduced gas temperature Trot versus the projected distance, d, to Feldt's star. Indeed, a 
good correlation between the temperature and the projected distance is revealed in Figure|l]a. 
In contrast, the correlation between Tj-ot and d becomes much more vague when the fiducial 
center is shifted about 2" in any direction away from Feldt's star. Together with the sufficient 
luminosity budget of the 05-type Feldt's star and the non-detection of other candidate Hii 
region powering stars within 2" from Feldt's star, we conclude that the majority of the energy 
for heating the dense gas enshrouding the UC Hii region G5.89 comes from Feldt's star. 

A power-law fit to the data indicates that the gas temperature scales as T^ot oc d~^''^^. 
Note that the deduced temperature Tj.ot{d) is actually a representative or averaged temper- 
ature for the gas along the line of sight. We suggest, however, that with a highly centrally- 
concentrated density distribution as shown in §5], the deduced temperature at small d is 
heavily weighted by high density, hence this represents the gas temperature at small r, 
where r is the distance to Feldt's star. Therefore the temperature profile T^ot{d) can be 
viewed as Ty-otii^)- The temperature profile of T oc r"*^'^*^ is shallower than the profile of T oc 



r measured toward hot-molecular core sources (e.g.. ICesaroni et al.lll998 



Beltran et al. 



20051 ). This disparity can be attributed to different dust grain properties and/or variant 
evolutionary stages. The temperature structure provides clues for the dust opacity index /?. 
Theoretical models found that, for example, the temperature structure of uniform dust shells 
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surroundinK a newly formed O star can be described by Tdust oc r~", where a ~ 2/(/5+4) 
( jWolfire fc Churchwelj 119941 ) . Modifications of tlie models with a centrally-concentrated 
density profile, however, are required to match the case of G5.89. 



4.3. Filling Factors & Column Density Profiles 

Figures [3]6 and [3]c show the derived filling factor and CH3CN column density, respec- 
tively. The CH3CN column density plotted in Figure [3]c, A'^tot,/=o.i, is the CH3CN column 
density scaled by a filling factor of 0.1 (see below). Figure shows the derived filling factor 
versus the projected distance to Feldt's star, and it clearly demonstrates that the molecular 
gas in the close vicinity of Feldt's star (i.e., d <6000 AU, mainly associated with N-lobe) has 
small filling factors ranging from 0.02 to 0.24, suggestive of a clumpy medium. The small 
filling factors can also be readily recognized from the low brightness temperatures (<20 K) 
of the hot (~100— 150 K) but optically thick low-iT components. For the gas further out, 
the filling factors become noticeably scattered with not only large values of 0.8—1 for the 
majority but also large uncertainties. Nevertheless, such large uncertainties in fact imply 
all the observed transitions to be optically thin. In such case, the filling factor, /, and the 
total column density, A^tot, become degenerate in the calculation and hence cannot be well 
determined. 

In short, the CH3CN gas toward G5.89 can be coherently considered as clumpy, with 
/ ~ 0.1.— 0.2. As shown in Figure Hie, we thus plot the CH3CN column density scaled by 
a filling factor of 0.1, A'^tot5/=o.i5 versus the projected distance to Feldt's star. Obviously 
■^tot,/=o.i decreases with increasing projected distance and the ratio between the inner {d 
^ 4000 AU) and outer {d ~ 18000 AU) region is more than a factor of 10. Note that if iVtot 
instead of N tot, f =0.1 is plotted, this ratio will be even larger because the inferred / is larger 
for gas further out. 



5. Density Structure 



The deduced CH3CN column density profile reflects the density structure of the molec- 
ular gas. For simplicity, we assume a spherical cloud core with a power-law density profile, 
n(r) ~ r~^, and a constant CH3CN fractional al 



^^^^^^ \ — 

Wilner et al. 


1994; 


Zhane et al. 


1998) 



We further adopt inner and outer boundaries at 4000 and 20000 AU, respectively, based 
on the size of the UC Hii region and the CH3CN emission. We then convolve the column 
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density with a two-dimensional Gaussian beam of 3.1" x 1.8", equivalent to our SMA beam. 
As shown in Figure Hie, we examine three model column density profiles calculated with 
power- law indices p of 0, 2, and 5.5. The adopted values for density at the inner edge are 
3.5x10^, 3.5x10^, and 3.5x10^ cm~^, respectively. Obviously, the modeled column density 
profile of p=0 (i.e., constant density) is too fiat to agree with the observational results (i.e., 
A''tot,/=o.i profile). Since the Ntot profile would be even steeper than A'"tot,/=o.i, the discrep- 
ancy between the profiles of observations and constant density can not be attributed to the 
adopted constant filling factor of 0.1. The observed Ntot,f=o.i profile can be well represented 
by the results modeled from a broken power law density distribution, with p = 5.5 for r 
< 8000 AU and p = 2.0 for 8000 AU < r < 20000 AU. Note that for the p = 2 profile, 
the gas density at r = 20000 AU is 1.4 x 10^ cm~^, more than sufficient to thermalize the 
observed CH3CN transitions. We therefore conclude that a centrally-concentrated density 
distribution is required to interpret the observational column density profile. 

Observations of massive star-forming cores indicated that the density power-law indices 
range from 1 to 2, while theoretical models also predicted similar power-law indices (e.g.. 



Mueller et al.ll2002l . and reference therein). The density profile of n ~ ^--ss jg significantly 
steeper than typical cases. Feedback processes such as the expansion of the Hii region 
may play a role to steepen the density profile. Since the observed CH3CN column density 
depends on the gas density and CH3CN abundance, ^(CHsCN), the steep CH3CN column 
density profile may also be a result of an enhancement of CH3CN gas due to the evaporation 
of grain mantles in the inner region, similar to the cases of abundance jump reported by 



van der Tak et all fl2000f ). 



Finally, we note that the temperatures as well as column densities are estimated simul- 
taneously under the assumption of one isothermal gas component along the line of sight. 
While the deduced temperature and density structures are in fact not uniform, they form 
a self-consistent set of solutions. That is, the highly centrally-concentrated CH3CN density 
profile supports the temperature profile, which in term validates the density profile itself. 
The derived "projected" profile therefore reflects the actual radial-dependent structure. More 
rigorous radiative transfer calculation that can account for all the involved parameters such 
as temperature, density, filling factors, and abundance, as well as observations with even 
higher angular resolutions will be desired to further refine the results. 
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for reading the manuscript. S.-Y. L., Y.-N. S., and H.-R. C. thank the National Science 
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Fig. 1. — The integrated K = 3 component of the CH3CN (12—11) emission shown in both 
red contours and color scales. Contour levels are —3, 3, 6, 10, 15, 20, 25, 30, 40, 50, 60, 
and 7 x 0.35 Jy beam ~^ km s^^. The gray contours represent the free-free emission at 
2 cm (iTang et al.ll2009l ). The dark ellipse at the bottom left denotes the synt hesized beam. 



The c rosses mark the positions of the sub-mm dust coi idensations repor ted by iHunter et al. 
( 120081 ) and the star marks the position of Feldt's star (jPeldt et al.ll2003l ). 
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Fig. 2. — The spectra of the CHs CN (12—11) hnes ta ken at the positions of the five sub-mm 
dust condensations identified by iHunter et al.l (120081 ) as well as Feldt's star. 
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Fig. 3. — The inferred rotational temperature (panel a), filling factor (panel b), and CH3CN 
column density (panel c) in color scales overlaid with the contours of the integrated CH3CN 
(12—11) K = 3 component. The CH3CN column density plotted in panel c, A'^tot,/=o.i, is 
the CH3CN column density scaled by a filling factor of 0.1. Caption as in Figure [H 



- 13 - 




2 4 6 8 10 12 14 16 18 

projected distance (kAU) 

Fig. 4. — The inferred rotational temperature (panel a), filling factor (panel b), and CH3CN 
column density (panel c) versus the projected distance, d, to the Feldt's star. In each panel, 
the overlaid greyscales represent the 1-a error in y-axis of the plotted data. In panel a, the 
best power-law fit. Trot oc is also plotted. In panel c, the blue, orange, and red crosses 

represent the modeled CH3CN column densities with the spherical density distributions of 
power-law indices, p, of 0, 2, and 5.5. See ^for the details. 



